HYDRODYNAMIC LIMIT FOR THE GINZBURG-LANDAU V0 
INTERFACE MODEL WITH A CONSERVATION LAW AND THE 
DIRICHLET BOUNDARY CONDITIONS 



TAKAO NISHIKAWA 

Abstract. Hydrodynamic limit for the Ginzburg-Landau V</> interface model with a 
conservation law was established in [5] under the periodic boundary conditions. This 
paper studies the same problem on the bounded domain imposing Dirichlet boundary 
conditions. A nonlinear partial equation of fourth order with boundary conditions is 
derived as the macroscopic equation, which is related to the Wulff shape derived by [2]- 



1. Introduction 

The Ginzburg-Landau V0 interface model determines stochastic dynamics for a dis- 
cretized hypersurface separating two microscopic phases embedded in the d+1 dimensional 
space. The position of the hypersurface is described by height variables <fi = {<j)(x);x G T} 
measured from a fixed d- dimensional discrete hyperplane V. We will take V = Tn := 
(Z/iVZ) d when we consider the system on a discretized torus with periodic boundary con- 
dition, or T = Djs; C Z d when we consider the system on the domain with some boundary 
condition. Here, is a microscopic domain corresponding to a given macroscopic do- 
main D <Z M. d which is bounded and has a smooth boundary. We then admit an energy 
(Hamiltonian) for the interface <fi by 

H(<t>) = \ E V(cj ) (x)-0(y))+ VMx) - (fry)) 

\x-y\=l \ x ~y\=l 

with a potential V G C 2 (R). Note that we need to give a boundary condition {<f)(x); x G 
Z d \ -Djv} in order to define the Hamiltonian H. Once we introduce the Hamiltonian H, 
the dynamics of the interface can be introduced by means of the Langevin equation 

d(j) t (x) = -U x {(j) t ) dt + y/2dw t (x), x G T, (1.1) 

where {wt(x);x G T} is a family of independent copies of the one dimensional standard 
Browninan motion, and U x (<p) is defined by 

\x-y\=l 



The hydrodynamic scaling limit has been established for the dynamics governed by (11. ip 
with T = Ttv in [5], for the dynamics with T = Dn and the Dirichlet boundary condition in 
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[7] . In both cases, the macroscopic motion is described by the nonlinear partial differential 
equation 

J^(*,0) = div{(V<r)(V/i(t,0)) 

d 



with an appropriate boundary condition. 

The dynamics (II. 11) can be regarded as the model corresponding to the Glauber dynam- 
ics in the particles' systems. Let us introduce the model corresponding to the Kawasaki 
dynamics in the particles' systems as follows: 

d(j) t {x) = -(-A r )[/.(0 t ) dt + y/2dw t {x), xET, (1.2) 

where Ar is the discrete Laplacian on T defined by 

yeT,\x-y\=l 

and {wt(x);x G T} is a family of Gaussian processes with mean zero and covariance 
structure 

E[w t (x)w t (y)] = -A v (x, y)t As, x, y e T, t, s > 0. 

We note that the dynamics f ll.2p preserves the sum Ylxtr ^tix), which can be regarded 
as the volume of the phase under the interface. 

The main purpose in this paper is to establish the hydrodynamic scaling limit of 4>t 
determined by ( 1 1 . 2 p under the Dirichlet boundary condition 

(j> t {x) = 0, x e Z d \ D N , t > 0, 

and to clearfy the relationship between the macroscopic motion and "Wulff shape" studied 
by [2]. The main result in this paper is that, under the scaling iV 4 for time while iV for 
space, the macroscopic motion corresponding to <fit is described by the nonlinear partial 
differential equation with Dirichlet boundary condition 

d -h(t,6) = -Adiv|(Vo-)(V/i(t,0)) 



dt 



i=l j=l J ' 



;i.3) 



h{t,9) = 0, 6eD c ,t>0. 



We note that the equation ( 11. 3p is the gradient flow in an affine space in (if 1 (/}))* with 
respect to the energy functional 



E D (h) = I a{Vh{6))dO 
Jd 
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with the Dirichlet boundary condition h\ D c = 0. Here, a : M. d — > M is a function called 
"surface tension," which gives the local energy of macroscopic interface with tilt u £ M. d , 
see [m for precise definition. The functional is called "total surface tension," which 
gives the total energy of interface h. We note that total surface tension T,d is the rate 
functional for the large deviation principle under the static situation, see [2]. Taking 
T = Ttv instead of Dn, the large scale hydrodynamic behavior has been studied by [6j. 

We should also mention the relationship between the equation (II. 3p and the Wulff 
shape discussed in [2] . As an application of the large deviation principle, the macroscopic 
height variable h N under the equilibrium state (Gibbs measure) conditioned on the total 
volume converges to the macroscopic interface so-called "Wullf shape," the solution of the 
variational problem 



I. J d ) 

as N — > oo, where v is the limit of the volume rescaled by N~ d . We emphasize that the 
solution h(t) for (II. 3p converges as t — > oo, and the limit coinsides with the solution of 
the variatinal problem (II. 4p . Indeed, the macroscopic motion described by (II. 3p relaxes 
the total energy E^> and attains to the Wulff shape as the limit t — > oo. 

Before closing the introduction, let us give briefly the organization of this paper. In Sec- 
tion [2J we formulate our problem more precisely and state the main result. In Section [3J, 
we study several properties of the macroscopic equation (11.31) and its spatial discretization. 
In Section HI we show that a translation-invariant stationary measure for the dynamics <f) t 
on the infinite lattice need to be canonical Gibbs measure corresponding to the Hamilton- 
ian H. Combining the above with the known result in [B], we have the characterization 
of the family of translation-invariant stationary measures. In Section we derive the 
macroscopic equation (II. 3p from the stochastic dynamics (II. 2p . after establishing several 
estimates for stochastic dynamics (12. ip . 



2.1. Model. Let D be a bounded, connected domain in M. d with a Lipschitz boundary. 
For convenience, let D contain the origin of M. d . In order to approximate profiles with 
the Neumann boundary condition by discrete ones, we will only consider D satisfying the 
following: 

Assumption 2.1. Let Dn = ND PI Z d . We assume that there exists a constant C > 
independent of N such that 

df) N (x,y) < C 2 , x,y E D N , \x - z\ < 2, \y - z\ < 2 
for every N >1 large enough and z £ Z d \ Dn, where is the ordinary graph distance 



Remark 2.1. Assumption 12. II is satisfied if the domain D is convex. 

Let us introduce the discretized microscopic domain corresponding to D. To keep 
notation simple, we shall consider Dn C Z d defined by 




(1.4) 



2. Model and main results 



on Dn- 



D n = {x £ D N ; B(x/N,5/N) c D} 
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where B(a,l) stands for the hypercube in M. d with center a and side length /, that is, 

d 

B{a,l) = 1 [[[a i -l/2,a i + l/2). 

i=i 

On Dn we consider the dynamics governed by the following stochastic differential equa- 
tions (SDEs) 

d<f) t {x) = -(-A DN )U.((j>t)(x) dt + \f2dw t (x), x G D N , (2.1) 
with the Dirichlet boundary condition 

(j) t {x) = 0, x eZ d \D N (2.2) 
and initial data 4>q, where U x {4>) in the drift term is defined by 

um--=-^M= E v'M*)-M) ( 2 - 3 ) 

''- \x—y\=l 



for (p G R Djv and x G Dn- Here, A^ is the Laplace operator on D^, that is, A^ is 
define by 

a Dn ^(x) = &(y) - M*)) 

2/G-Djv; \x-y\=l 

for ip N : R Djv — > R. The process w)t = {i2; t (x); a; G -D^} is a family of Gaussian processes 
with mean zero and covariance structure 

E[w t (x)w s (x)} = (-A DN )(x,y)t As, x,y G D N , t,s>0. 

Note that a stochastic processes Wt satisfying the above can be constructed by 



W t (x) = y/-A DjT W t (x), 

where ^—A Dn is the square root of — A Dn and {w t {x)\ x G D N } are a family of indepen- 
dent one dimensional Brownian motions. For convinience, we extend Wt to the process on 
"L d by putting Wt(x) = when x G Z d \ D^. 

Through out this paper, we always assume the following condition on V: 

Assumption 2.2. The function V : R — > R satisfies the conditions as follows: 

(1) V G C 2 (R). 

(2) V is symmetric, that is, V(rj) = V(—rj) holds for all r\ G R. 

(3) There exist constants c + , c_ > such that 

c- < V"{ri) < c+, ?7 G R 

holds. 

We regard (12. ip as the model describing the motion of microscopic interfaces and in- 
troduce the macroscopic height variable h N by scaling N A for time while N for space: 

h N (t,9) = £ N- l <j> NH (x)l B{x/N , l/N) (e), 9 G R d , 

x& d 

where <p t = {0t(x); x G Z d } being the solution of §T[§ with lj23]> . Note that the suitable 
scaling is not the diffusive one. 
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2.2. Main Result. The main result in this paper is the following: 

Theorem 2.1. We assume Assumptions \2.1\ and \2.2\ We furthermore assume that the 
sequence of initial data 0o — f or P-ip satisfies 

lim £11^(0) - h \\ 2 H UD)* = ( 2 - 4 ) 

with some h £ L 2 (D), where h N (0) is the macroscopic height variable corresponding to 
4>q . Then, for every t > 0, h N (t) converges as N — >■ oo to h{t), which is the unique weak 
solution of the partial differential equation (11.31) with initial data h$. More precisely, for 
every t > 0, 

lim E\\h N (t)-h(t)\\ 2 Hl{Dr =0 (2.5) 

holds. 



3. The macroscopic equation and its discretization 

In this section, we shall focus our attention on the limit equation (11.31) and its discretized 
version. The arguments in this section highly depend on the properties of the surface 
tension a established in [2] and [5]. 

3.1. Sobolev space. As we will see, our computation is based on i/~ 1 -norm. Before 
starting calculations, we introduce if~ 1 -norm and its discretization. Since the solution of 
(ll.3p should satisfy 

f h(t,6)d8= [ h (9)d6, t>0, 

J D J D 

which means that the actual state space for (II. 3p is not linear. For convinience, let us 
mainly consider the time evolution on the tangential space. We denote the dual of the 
Sololev space H l (D) by if 1 (£>)*. The equation (11.31) will be actually solved at 

H = {he H\D)*; (M) = 0}, 
where (•, •) denotes the duality pairing between H 1 (D)* and H 1 (D) such that 

(f,g) = (f,g)^(D) 

for every / £ L 2 (D), where (-, -)l 2 (d) denotes the inner product in L 2 (D). 

We shall next introduce the Green operator G for the Laplacian with zero Neumann 
boundary condition. For h £ H, we denote the unique solution g £ if x (Z)) of the elliptic 
equation 

(Vg,VJ) = (h,J), JeH\D), 
= 

by Gf. Let us define the bilinear form (•, on H by 

Oi, h 2 ) H = (VGhx, VG7i 2 ), h h h 2 £ H. 

We then have that H is the Hilbert space with the inner product (•, -)h, and that the 
norm associated to (•, •)# is equivalent to the ordinary H 1 (D)* norm restricted to H . 
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Next, let us define the discrete version of the Hilbert space H define above, analogously. 
For the step function / : M. d — > M. with mesh size 1/N, that is, / has the following 
representation 

f N (6) = J2N- 1 i> N (x)l B(x/N , 1/N] (9) 

x&L d 

with i/j n G M z< \ we define \\f N \\^ N by 

\\f N \\\ N = Yl (^(*) - (^x-A^r 1 ^*) - <^» + N-™-*{r)\ 

where 

x£D N 

Note that the inverse of the Laplacian (— A^) -1 can be defined as the linear operator 
from 



N . 



E *(*) = °} 

t£D n ) 



x£D N 

to itself, by regarding A Dn as 1:1 onto map from stf N to ^/ N . We can easily see that there 
exists a constant C > independent of N such that 

imi^W < CII/^ll-Liv- (3.1) 
holds for every step function / with mesh size 1/N satisfying / | D o = 0. We will 
frequently use (13.11) to establish a priori bounds and so on. 

3.2. Precise formulation for macroscopic equation. We shall give the precise mean- 
ing of the solution of ( II. 3p . Let us introduce a triple V C H = H* C V* by H introduced 
in Section I3TT] 



V= j/i G Hl{D); j h(t,0)d9 = 0^ 



and the dual space V*. Note that the space V introduced above is closed subspace of 
Hq(D) and therefore V is reflexive. We denote the duality relation between V* and V by 
v*(-, -)vi which satisfies 

v*{f,g)v = U,9)h 

for / G H and g G V. Let us consider the nonlinear fourth order differential operator 
Af.V^V* defined by 

A f (h) := -Adiv((Vtr)(V/i + V/)) , ft G V 

for given / G Hq(D). We recall that the surface tension a : M d — >• 1R is C 1 -class and it 
satisfies 

c_\u-v\ 2 < (Va(u) - Va(v)) ■ (u - v) < c + \u-v\ 2 , u,v eR d . (3.2) 

Definition 3.1. A function h = h(t,9) is called the solution of (11. 3p with initial data 
h G H X {D)* if there exists a function / G Hq(D) with 

(ho,l) = (f,l) (3.3) 
such that the function hf := h — f satisfies the following conditions: 
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(1) hf : [0,T] — > V* is absolutely continuous and 

h f G C([0, T]; H) n L 2 (0, T; 1/) n W l ' 2 {[d, T], V*). 

(2) ^(0) = /i - /■ 

(3) /i/ satisfies 

±h f = A f (h f (t)) 
in for almost every t G [0, T]. 
Remark 3.1. The third condition is equivalent to 

h f (t) = hf(0)+ [ A f (h f (s))ds 
Jo 

in V* for almost every t G [0, T]. Roughly saying, the above is equivalent to 

= h{0) + / A(/i(s))ds, 

where A is the nonlinear fourth order differential operator defined by 

A(h) := -Adiv((Va)(V/i)). 

The first aim in this section is to see the existence and uniqueness for our equation. 

Theorem 3.2. For every h G H l (D)* , there exists a unique solution of (11.31) . 

In order to apply a general theory of nonlinear partial differential equations, we prepare 
the following lemma: 

Lemma 3.3. For every f G Hq(D), the nonlinear operator Af : V — > V* satisfies the 
following: 

(1) A is monotone, that is, 

v * (Ajih) - A f (h 2 ), h x -h 2 ) v < 0, h u h 2 e V. 

holds. 

(2) A is demicontinuous , that is, Aj(h) is continuous under the weak topology ofV*. 

(3) There exist constants C\,C 2 ,C^ > such that 

v*(Af(h),h) v <-C 1 \\h\\ 2 v + C 2 , heV, (3.4) 

\\A f {h)\\ v .<C 3 (\\h\\v + l), heV (3.5) 

holds. 

Proof. We at first note 

v* (A f (h),g) v = -((W)(V/i + V/), Vg) L 2 {D) (3.6) 

for h, g G V, and 

v^Afih) - A f (h 2 ),g) v = -((W)(V/n + V/) - (W)(V/i 2 + Vf),Vg) L 2 {D) (3.7) 

for hi, h 2 , g G V. It is easy to see (1), by applying (I3.7p with g = hi — h 2 and using the 
convexity of a, see (13.21) . We can also obtain (2), since we have 

\v*{A f (hi) - A f {h 2 ),g) v \ < c+HV/ii - Vh 2 \\ L 2 {D) \\Vg\\ L 2 {Dh h u h 2 ,ge V (3.8) 
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from (13. 2p and (13.71) again. Moreover, we can obtain (13. 5p of (3) also, because the rela- 
tionship ( 13. 6 p implies 

\\A f (h)\\y < \\(Va)(Vh + Vf)\\ L 2 {D) <c + \\h\\ v + c + \\Vf\\ LHD) 

for heV. 

As the final step, we shall show (13.41) of (3). Using (13.61) . we obtain 
v . (A f (h), h) v = -((Va)(Vh + V/), Vh) L2(D) 

= -((Va)(Vh + V/), Vh + Vf) LHD) + ((W)(V/i + Vf), Vf) LHD) 
< -c_||V/i + V/||| 2(D) + c+\\Vh + Vf\\ L 2 (D) \\Vf\\ L 2 (D) . 



Here, noting 



\Vh + Vf\\h iD ) > \\\Vh\\h {D) - ||V/||i 2(D) , 



and 



c+\\Vh + V/|| L 2(d) 



nv/ir 



L 2 (D) 



we obtain 



v*(A f (h), h )v < -\c-\Wh\\ L * {D) + + c + + c_) || Vf\\ 2 L2{D) . 
Applying the Poincare inequality for h G V C Hq(D), we conclude the desired bound 

dS3D- □ 

Proof of Theorem \3.2 . Using Theorem 4.10 of p] and Lemma [3.3 1 for every initial data 
h G and every auxiliary function / G Hq(D) satisfying (13.31) . we obtain the 

existence and uniqueness of hf satisfying conditions (l)-(3) in Definition 13. 1[ Especially, 
this shows the existence of the solution of (11.31) in the sense of Definition 13.11 

Let us show the uniqueness of of the solution of (II .3p in the sense of Definition 13.11 
Take a two solution W- 1 ' and with a common initial data ho, and let f\ and f^ be 
auxiliary fuctions associated to and h^ 2 \ respectively. Noting 

A fl (h)=A f2 (h-f 2 + fi), heV, 

— fi = ft"? + /2 — fi satisfies conditions (l)-(3) in Definition 13.11 with f\. The 
uniqueness of hf for given / implies h^ — fx — — fi, which shows h^ = h^ 2 \ □ 

3.3. Regularization for the macroscopic equation. Let us introduce the regulariza- 
tion of a and the corresponding partial differential equation, which plays key role in the 
proof of Theorem 12. 1[ Note that such regularization is not needed once one can solve 
C 2 -regularity of a. 

Let p G C^(M. d ) be non-negative, symmetric, suppu C {u G M. d ; \u\ < 1} and 
f Rd p{x) dx = 1. For < S < 1, we define ps by 

P5 = 5- d p(5~ 1 u), ueM. d 
and the regularized surface tension a 5 by the mollification of a: 

a 5 {u) = a * p 5 (u), u G M. d . 
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Note that the regularized surface tension a 5 again satisfies the bound (13. 2p . that is, 

c„\u-v\ 2 < (Va s (u) - Va 5 {v)) ■ (u - v) < c + \u-v\ 2 , u,veR d (3.9) 
holds. Moreover, since Va is continuous, a 5 approximate a in the following sense: 

\X7a s (u) - Va(u))\ < c + S, uER d . (3.10) 
Using o~ s defined above, let us consider the nonlinear fourth order differential equation 



^h 5 (t,9) = -Adiv{(V</)(V/i(t,0))}, 9 E D, t > 
h s (t,6) = 0, 6eD c ,t>0. 



The equation (13.111) can be formulated by similar way to Definition 13.11 Since we have 
(13. 9p . Theorem 13.21 can be applied to (13. lip also. We therefore obtain that the equation 
(13. lip has a unique solution. Furthermore, we get the following proposition, which implies 
that the solution of (ll.3p can be approximated by the solution of (13. lip . 

Proposition 3.4. Let h and h s be the solutions of (ll.3p and (13.111) . respectively. If the 
initial datum are common, we then have 

\im\\h(t)-h 6 (t)\\ HHDr = 0. 
Proof. Let us fix an auxiliary function / e Hq(D). By Definition 13.11 we obtain 

\\h\t) - f\\ 2 H = \\h\0) - f\\ 2 H - 2 [\va 5 (Vh 5 (s)) - Va 5 (Vf), Vh 5 (s) - Vf) L 2 {D) ds 

Jo 

-2 / (Va\Vf),Vh\s)-Vf) L 2 {D) ds 
Jo 

< \\h S (0) -f\\ 2 H - 2c_ / \\Vh s (s) - Vf\\ 2 L2{D) ds 







+ 2 / \\Va d (Vf)\\ L2{D) \\Vh d (s)-Vf\\ L 2 {D) ds 
Jo 

< \\h s (0)-f\\%-c- [ \\Vh s (s)-Vf\ 

Jo 



I 2 ds 



+ CZ 1 [ || Va s (Vf) \\ 2 L , {D) ds 
Jo 

from (13. 6p and (13.91) . This shows 

sup/ \\h s (s)\\ 2 HHD) ds < oo. (3.12) 

5>0 Jo 
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By similar calculation to the above, we obtain 
\\h(t)-h'(t)\\ 2 B 

= -2 [ (V(t(V/j(s) + V/) - Vff(V^(s) + V/), V/j(s) - V^(s)) L 2 (fl) (is 
Jo 

-2 / (Va(Vh 5 (s) + V/) - Va 5 (Vh s (s) + Vf),Vh{s) - Vh 5 (s)) L2{D) ds 
Jo 

< 2 I c + 5 [ \Vh(s,8) - Vh s (s,8)\d8ds. 
Jo Jd 

Combining above with (I3.12p . we get the conclusion. □ 

We can show the following proposition by similar argument to the proof of Proposi- 
tion E31 

Proposition 3.5. Let both of h and h be the solutions of (11.31) . We then have 

\\h(t)-h(t)\\ 2 H <\\Ko)-h(o)\\H, t>o. 

3.4. The discretization for the macroscopic equation. In order to introduce the 
discretized equation corresponding to the regularized macroscopic equation (13. lip , let us 
introduce several notations. We define the finite difference operators by 

Vif(P) = N(f(e + ei/N)-f(0)) t 

V?'*f{0) = -N{f{0)-f(p-e i /N)), 
V"/(0) = (< /(*),..., V?/(0)), 

d 

dw N g(6) = -J2^9m 

i=l 

for / : R d E, g = (g l )i< i < d : R d -> R d , 1 < i < d and 9 e R d , where e { G Z d is the i-th 
unit vector given by (ei)j = 5ij for 1 < % < d. We also define the discretized Laplacian 
with the Neumann boundary condition by 

d 
i=l 

for / : M. d — > R, where the domain is defined by 

D N = |J B(x/N,l/N). 

Note that indicator functions appearing above are corresponding to the range of x's where 
the sum (12.31) is taken. With these notations the discretized PDE for (13.111) reads a system 
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of ordinary differential equations 

' ^h N,5 (t, 6) — —A N k N ' 6 

:= dw N {(Va s )(V N h N ' 5 (t))}(6), 9 G D N , 
h N ' s (t,6) = 0, 9£D N . 

We recall that a 5 G L7°°(M d ) and it satisfies 1ET5]) . The equation fl3~T3]) will be solved with 
the initial data given by 

ho' s {6) = N d [ h (9) dO, 9 e B{x/N, 1/N), x G Z d . (3.14) 

JB(x/N,l/N) 

where ho G L 2 (M. d ). Since the initial data is a step function, the solution h N (t,6) is 
also a step function, that is, 

h N > s (t, 9) = h N ' 5 {t, x/N), 6 G B(x/N, 1/N), x G Z d . 

Though the equation (13. 13[) relies on the value of k N ' s in only, let us extend k N ' s to 
the domain large enough for convinience. We define diD^ by 

diD^ = {x G Z d \ Dm] dist Z d(x, Dn) = i}, 

where dist Z d is the graph distance on Z d . We define the value of on dkD^ recursively, 
by 

k N >\x) = (#{yeD N ;\x-y\ = l})- 1 ]T k N > 5 (y), x G 8 1 D N , 

y€D N ;\x-y\=l 

and 

k N > s (x) = (#{yed l „ 1 D N ;\x-y\ = l}r 1 ^ k N > 5 (y), x^d l D N . 

y£d i - 1 D N ;\x-y\ = l 

3.5. A priori bound for the discretized equation. In this section, we establish a 
priori bound of the solution h N ' s of (13. 13p . which is uniform in the mesh size N. To do 
so, we introduce an auxiliary function similarly to Definition 13. II Let us take the function 
g G C^°(D) which satisfies the following: 

(1) g{x) > holds for every x G D. 

(2) J D g(B)d0 = l. 

We can then take N > 1 enough large such that suppg C D N for every N > N . Let us 
establish a priori bound for h N with N > N . 
Using g introduced above, we define ( N by 

( N ( x ) = N d+1 [ g(9)d8, x e Z d 

JB(x/N,l/N) 

and g N by 

g N {6) = N~^ N {x)l B{x/N , l/N) {e), 6 G W l . 

x£Z d 
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We then have the following bound: 
sup sup sup { \g N (x/N) | -f 

N>ll<i,j<d xe D^ 

where c g is the constant defined by 



sup sup sup {\g N (x/N)\ + \V?g N (x/N)\ + |Vf Vf g N (x/N)\} < c 

N>ll<i,j<d xe D^ 



(3.15) 





dg 




d 2 g 




09, 


+ 

oo 


d9 i 9 j 


oo y 



(3.16) 
(3.17) 



l<i,j<d 

Using ( N and g N introduced above, we define ip N and f N by 

^j n ( x ) = v( N (x), xez d 

and 

f N (x)=vg N (9), 9eR d 
respectively, where v is the "volume" of ho, that is, 

h^\9)d9= (ho,l) 

D 

for the initial datum h$ for (I3.13p . Note that the right hand side of the above does not 
depend on the choice of S. We then have that sequences {ip N } and {f N } satisfy following 
properties: 

(1) For every N and x G Z d \ D N , ip N (x) = f N (x/N) = holds. 

(2) For every N, the relationship 

4> N (x) = N j f N {9)d9 

holds. 

(3) The following bounds hold: 

sup sup sup{\f N {x/N)\ + \V?f N {x/N)\ + \V?Vff N {x/N)\}<c 9 

N>1 l<i,j<d x(E D^ 

Proposition 3.6. There exist constants C\,C2 independent of N,5 such that 

11^)11-!,* + N~ d f \\V N h N > 5 {s)\\% {D) ds < C^m-^N + C 2 (l + t) 



(3.18) 



holds for every t > 0. 
Proof Differentiating \\h N > 5 (t) - f N \\ 2 
d „ T 



1 N in t, we obtain 



dt 



\h N '\t)-r 



-1,7V 



= -2N- d V N h N < s (t,x/N)-Vcr 5 (V N h N ' 5 (t,x/N)) 
+ 2N~ d V N f N (x/N)-Va s (V N h N > s (t,x/N)) 
< -2c_N- d J2 \^ N h N ' S (t,x/N)\ 2 

x£D N 

+ 2N~ d y N f N (x/N) ■ Va 5 (V N h N ' s (t, x/N)), (3.19) 
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where Dm is defined by 

Dn = {x £ Z, d ; there exists y £ such that \x — y\ < 1}. 
Here, we have used the summation-by-parts formula 

oc(x/N) div N P{x/N) = - ^ V N a(x/N) ■ 0{x/N), (3.20) 

where a : M. d — > M and (3 = (Pi)f =1 : M d — > M. d are arbitrary functions such that a 
and /3j (1 < i < d) are step functions with mesh size 1/N and a(x/N) = for every 
x £ Z d \ D N . 

The second term in the right hand side can be estimated in the following way: 



N - d J2 V N f N (x/N) ■ Va s (V N h N > 5 (t,x/N)) 

< -c__iV d \V N h N ' s (t,x/N)\ 2 + CN- 



D 



with a constant C > independent of N. We have used the properties of f N stated at 
the beginning of this subsection. Plugging the above into (13 . 19[) and integrating in t, we 
obtain 



\h N ' 5 (T)-f N \\\ N < ||^(0) -f N \\\ N 



\V N h N ' 5 (t)\\ 2 L2{D) dt + CN- d \v\ 



D 



N 



T. 



□ 



which implies the desired estimate, since H/^H-^jv is bounded unformly in N. 

We can improve the bound for 'V N h N,s if the initial datum is smooth enough. 
Proposition 3.7. We assume that 

h £ C™{D) (3.21) 
We then have the following unform bound: 

sup sup \\V N h N ' s (t)\\ 2 L 2, D) < oo 

N 0<t<T 

for every T > . 



Proof. Differentiating ^ x€ jj^ o~ s (V N h N ' s (t, x/N)) in t, we have 

a\V N h N '\t,x/N))= Vo-\V N h N > 5 (t,x/N))V N ^-h N > 5 (t,x/N) 



Ot 



x&Dp 



Y div N Va 5 (V N h N '\t, x/N))^-h N >\t } x/N) 
' ot 



k N > 5 (t,x/N)A N k N ' s (t,x/N). 



(3.22) 
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Here, we have used (3.8) in [7J, since 

dh N 

-Qj-(t,x/N) = 0, xeD c N 

holds. Since — A^v is non- negative definite, we obtain that the right hand side is non- 
positive. Dropping the right hand side and integrating in t, we have 

cr 6 (V N h N ' 5 (t,x/N)) < ° S (V N h N ' s (0,x/N)), 

xeD N x£D N 

which indicate the conclusion, since the function o satisfies (13. 2p . □ 

Let us establish the bound for k N in (I3.13p . in order to apply the argument in Section 3.3 
of 0. 

Proposition 3.8. We assume (I3.2ip . We then have the following bound: 



sup < sup 

N 0<t<T 



dh 



N,S 



dt 



-1JV JO 



T 



/A' 



dh 



N,8 



dt 



dt><oc. 



Proof. Noting 



d 2 h 



dt 2 



i,j=l < * 3 J 



and 



we have 



-z-(t,x/N) = 0, 



dt 



d_ 

dt 



dh 



N,S 



dt 



-1,7V 



d f)h N > 5 r) 2 rr 



i.i=l x&D N 
i=1 x£D N 



dt 



duiduj 



(V N h N ' s )\/f^^- 



TV 



dh 



7V,<5 



dt 



-(t,x/N) 



by performing the summation-by-parts several times. Integrating the both sides in t, we 
have 



dh 1 



dt 



2c_N~ 



-1,7V 



EE (v{ 

* =1 x£D N 



TV 



dh 



7V,<5 



ds 



-(s,x/N)) ds< 



dh 



TV 



dt 



"(0) 



-1,7V 



which implies the conclusion. 
Remark 3.2. By the definition of k N,s , we have 
dh N > s 



□ 



dt 



-1JV 



N~ d k N ' 5 (t,x/N)(-A N k N ' s )(t,x/N). 



xeDp 
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We therefore obtain 

sup sup N~ d k N ' s (t,x/N)(-A N k N ' 5 )(t,x/N) < oo. 



N 0<t<T 



by Proposition 13.81 

Remark 3.3. We need the smoothness of a 6 in order to obtain the uniform bound in N 
for 

N~ d Y, k N ' s (0,x/N)(-A N k N ' s )(0,x/N) : 

even if ho is smooth, for example, ho G C^°(D). We have C 1 -regularity of a and the 
Lipschitz continuity of Vcr, but such regularity is less than that we need. It is the reason 
why we consider the equation (13.13P with smooth a 5 instead of the original surface tension 
cr. 

As a direct consequence of Proposition 13.81 we have the following result. 

Corollary 3.9. We assume (13.211) . There exists a constant C > independent of N such 
that 

sup - h N \t 2 )\\\ N < C\t x - t 2 \ 

N 

holds for every < t±, ti < T. 

3.6. Uniform L 2 -bound for k N,s . In this subsection, we shall establish the uniform 
L 2 -bound for k N,s . Our goal in this section is the following: 

Proposition 3.10. Under the assumpotion (13.211) . we obtain 

sup sup \{k N ' 6 {t))\ < oo (3.23) 

N 0<t<T 

for every T > 0, where 

{k N - 5 (t)} — N~ d k N ' 5 (t jX /N). 

In order to show the above, we shall reduce our problem to that for solutions of elliptic 
equations whose mail term is linear. We at first note that the gradient of the surface 
tension is expressed by the expected value of V ! {rj{b))\ 

V l o-(u) = E^[V'(T 1 {(e^)% l<t<d, 

where is the i-th unit vector in M. d . For simplicity, we denote the directed bond (e^, 0) 
simply by again. We then have the following decomposition for Vcr: 

Va(w) = A(u)u + a(u), (3.24) 

where A(u) = (Aij(u)) is the d x d diagonal matrix defined by 

• i 



Aij(u) = 



V"(rj(ei) - Xm) d\ 



1 <i,j <d 



and a(u) = (a,i(u)) € M d is defined by 

Oi{u) = [V^) - Ui )} , 1 < i < d. 
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We remark that the matrix A(u) satisfies 

c_I < A{u) < c+I (3.25) 

uniformly in u, where I is the d x d identity matrix. Furthermore, we also remark that 
the vector a satisfies 

sup sup |aj(u)| < C a 

l<i<d U £M. d 

with some constant C a > 0, which is a simple consequence of Brascamp-Lieb ineqality, 
see [2] and [5]. By (I3.24p and the definition of a s , we also obtain the decomposion for 
Va 5 as follows: 

Va\u) = A\u)u + a\u), 

where 

A$ j (u)=A ij *p s (u), l<i,j<d 

and 

af(u) = dij * ps(u) + / Aaiu — v)vips(v) dv, 1 < i < d. 

We note that A s is diagonal and satisfies (13.251) again, and that a 5 satisfies 

sup sup sup | of (it) | < C' a 

5>0 l<i<dui 



with some constant C' a > 0. Putting 

A N ' s {t,B) = A s (y N h N > 5 (t,6)), 
a N >\t,6) = a\V N h N >\t,6)), 

WC SCO that h N > s (t) satisfies 

k N ' 5 (t) = div N (A N > s (t)V N h N > 5 (t))) + dw N a N > 5 (t), t > 0. (3.26) 

Let us regard (13.261) as the elliptic equation for given A N ' s (t), a N ' s (t) and k N ' s , which has 
a unique solution. Since the main term of (I3.26P is linear, we can write h N,s as the sum 
of profiles. For t > 0, let hi' 5 (t) be the unique solution of 

~k N ' 5 (t) = dw N (A N > 5 {t)V N h^ 5 {t))) + divjv a N ' 5 (t) (3.27) 

with the Dirichlet boundary condition 

h" > 5 (t, x/N) =0, i£Z d \ D N , 

where u N)5 {t) is defined by 

~k N >\t) := k N ' 5 (t) - (k N ' 5 (t)}. 

Furthermore, we let ^^^(t) be the unique solution of 

(k N ' s (t)) = dw N (A N ' S (t)V N h^' 5 (t))) (3.28) 

with the Dirichlet boundary condition similarly to . Note that the original h (t) can 
be recovered by 

h N (t) = h?(t) + h%(t). (3.29) 
We shall at first show the bound for h± . 
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Proposition 3.11. There exist constants C\,Gi > independent of N such that 

liv^wili^ + \M (t)\\l HD) < c.fk^mi^ + c, 

holds for every t > 0. 

Proof. Multiplying the both side of (I3.27P by hi(t) and taking the sum over D^, we have 
J2 k N ' 5 (t,x/N)h?(t,x/N) = h?(t,x/N)div N (A N (t)V N h?(t)))(x/N) 

+ h?(t,x/N)div N a N (t)(x/N). 

Dividing the both side by N d and performing the summation-by-parts, we obtain 
N~ d ~k N ' 5 (t,x/N)h?(t,x/N)=N~ d J2 V N h?(t,x/N) ■ A N (t)V N h?(t,x/N) 

+ N~ d Vh?(t jX /N)a N (t,x/N) 
xeD N 



<-lc- \V N h?(t,x/N)\ 2 + ^N d \D N \ 



2 

xG-Djv 

We have used (I3.25P and h\{t,x/N) = for x G Z d \ D N . Using the Poincare inequality 

\\hnml><C\\V N h»(t)\\% 
for a constant C > independent in N, we have 

ic_||V^f(t)||| 2 < 2 7 ||fc^(t)||£ a + 2^\\hnt)\\h + ^N*\D N \ 

< 2 7 ||^(t)||i 2 + 2 7 - 1 L7||V/^(t)||! 2 + ^^1^1 
for every 7 > 0. Choosing 7 = 4C/c_, we conclude 

\c^\\v»h»{t)\\b < ^fk N ' s (t)\\h + ^n*\d n \. 

Applying the Poincare inequality to the above, we also obtain the bound for 11^(^)11^2- D 

Since we now have the nice bounds for h± (t) and h (t), we have one for h^ (t) also. 
We shall show the following proposition which says that the bound for h^ (t) implies the 
bound for (k N ' s (t)). 

Proposition 3.12. For a E R, let h^ai^) be the solution of 

a = div N (A N (t)X7 N h%Jt))) (3.30) 
with the Dirichlet boundary condition. We then have 

a? < C\\V N hZJl 2{D) . (3.31) 
with a constant C > independent of N,t and a. 
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Proof. Multiplying (I3.3(jp by h^ a and taking sum over D^, we have 

aN' d = N ~ d E h% a (t,x/N) div N (A N (t)V N h» a (t,x/N)) . 

Performing the summation-by-parts at the right hand side, we obtain 

aN~ d x/N) = -N~ d V%(i, x/N) • A*(t)V%(t, nr/JV). (3.32) 

We have used h^it) satisfies the Dirichlet boundary condition. Noting that h^ a can be 
expressed by 

h 2,a = ah 2,l 

because the right hand side of (13.301) is linear, we get 



a 



l N~ d £ h^(t,x/N) 



= N~ d Y V%(t,i/JV) • A N (t)V N h» ttt (t,x/N) 

x£D N 

<c„N~ d J2 \V N hUt,x/N)\ 2 



x£Dm 



from fl3T25l) . 

Here, once we have 



inf 

N 



> C 



(3.33) 



N~ d £ hlS, x IN) 

with a constant c > independent of iV and i, we immediately obtain the conclusion. We 
shall therefore show (I3.33p . Using f !3.32j) and (13. 25 p . we get 



> r \\V N h N II 2 



(3.34) 



For the function g introduced at the beginning of Section 13. 5[ we obtain 



1 = N~ d Y 9 N (x/N) 

x£D N 

= N~ d Y V N 9 N (x/N) ■ A N (t)V N hl 1 (t,x/N) 



< civVIlWIvXiWIUw 



with a constant C > by using (I3.25p . Combining the above with (13 . 1 5[) and (I3.34p . we 

get 



N~ d E hUt,x/N) 



xeDp 



> 



C' 2 c 2 g ' 



which shows (13. 33 p . 



□ 
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Once we obtain Proposition 13.111 and Proposition 13.121 we can easily show Proposi- 
tion GD3B Noting 

\\V N h»(t)\\bw < 2|| W(*)||i» (B) + 2||V^ (t)\\l HD) , 

we have 

\\<nNi.N/j.\U 

2(D) 



SUp SUp || V JV /i^(t)||_L2( D ) < OO 



N 0<t<T 

by using Proposition 13.71 and Proposition 13.111 On the other hand, since we have 

\(k N ' s (t))\ 2 <C\\V N h»(t)\\ L * {D) 
from Proposition 13.121 we obtain (I3.23P and therefore Proposition 13.101 
Corollary 3.13. Under the assumpotion (13.211) . we have the following bound: 

sup sup {\\k N > 5 {t)\\ LHD) + || V N k N > 6 (t)\\ L 2 (D) ) < oo. (3.35) 

N 0<t<T 

Proof. From Assmption 12.11 and the definition of k N ' s as in Section 13.41 we can easily see 

P^'(t)IU»( D )<Ci J2 \k N '\t,x/N)\ 2 

xGD N 

and 

II V N k N ' s (t) \\l HD) <C 2 N~ d k N > s (t,x/N)(-A N k N > 5 )(t,x/N) 

x£Dpf 

with constants C\,Ci > independent of N. These inequalities and Proposition 13.101 
imply fl3~35D . □ 

Once we have Corollary 13.131 we can also obtain uniform L p -bound for W N h N . This 
uniform bound plays the key role in the derivation of PDE (11.31) from the height variable 
h N . 



Proposition 3.14. We assume (13.211) . We then have the following bounds: 

p 

Lp(D) 



sup sup \\k N ' s (t)\\ p LP(D) < oo 



N 0<t<T 

and 

7NZN,5u\\\P 



sup sup ||V w /i JV ' d (t)||^ (D) < OO 

N 0<t<T 



for some p > 2. 



Proof. Combining Corollary 13. 131 with the similar argument to the proof of Proposition 1.4 
in [5], we obtain the first assertion. Applying the argument in Section 3.3 of [7] to (13.271) . 
we also obtain the second assertion. □ 

As an application of Proposition 13.141 we can obtain the following identity correspond- 
ing to the oscillation inequality in [2]. This also plays the key role in the derivation of 
PDE (11.31) from the height variable h . 
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Proposition 3.15. We assume f)3.2ip . For the solution h N,s of (I3.13P and e6Z d such 
that \e\ = 1, we have 

T 



(3.36) 



lim N~ d [ V \V N h N ' s {t,x/N + e/N) -V N h N ' s {t,x/N)\ 2 dt = 

Proof. Take 1 < % < d arbitrary. From (13.1 3p . we can split 

F N (t) :=N^ 2 V?k N < s (t,x/N)V?h N ' 5 (t,x/N) 

x£Dm 

into three terms as follows: 

d 

F ^ t)=N -d-2 Y,V?Va(V N h N >\t,x/N))V?Vfh N '\t,x/N) 
xeTT^ ■ ?=1 

_ N -d-2 Vf diw N Va(V N h N ' 6 (t,x/N))Vfh N ' 6 (t,x/N) 

x£D^n(D N - ei ) c 

+ N~ d ~ 2 Yl V?k N ' 6 (t, x/N)V?h N ' 5 (t, x/N) 

x£D^n(D N - ei ) c 

Here, we have used 

V N h N ' 5 (t, x/N) = 0, xeZ d \TJ^. 
Using (GOD, we have for Ff(t) 

F?(t) < -c_N~ d \V N h N ' s (t,x/N + ei/N) - V N h N ' s (t,x/N)\ 2 , 
xeDpf 

which is nothing but our target. From now on, we shall show the remaining terms vanish 
when N — > oo. 

For F^(t), since we have 

rT r T 



L 2 (D) at ' 



\F N (t)\<2N~ 2 f \\V N k N > 5 (t)\\ 2 L2(D) dt + 2N- 2 [ \\V N h 
Jo Jo 

by Schwarz's inequality, we obtain 

lim [ T \F N (t)\dt = 

from Proposition 13. 71 

We shall next establish the bound for and F^ . To do so, we shall at first make an L? 
bound for \7 N h N on Bn = \ D^. Choosing r, q > 1 such that r = p/2, l/r+ 1/q — 1 
and applying Holder's inequality, we obtain 

N~ d \V N h N '\t,x/N)\ 2 < ~(N- d \B^\J /2q N- d \V N h N ' s (t,x/N)\ p 



N 



1 / \ 1/2 

-(N- d \B N \) . (3.37) 
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For F 2 (t), since we have 

\F 2 N (t)\<N~ d \V N h N ' s (t,x/N)\ 2 , 

we obtain 

lim / \F 2 N (t)\dt = 
from Proposition 13.141 Also for F^(t), we have 

\Ff(t)\ < 2a(N)N~ d \V N h N (t,x/N)\ 2 + 2a(Ny l N~ d | W'(t, x/N) f 

for an arbitrary sequence {a(N)} of positive numbers. Choosing a(N) = N e with e > 
small enough, we obtain 

lim [ T \F 3 N (t)\dt = 0. 

from Proposition 13.141 Summarizing above, we conclude ( I3.36p . □ 

3.7. Convergence of the solution for the discretized equation. In this subsection, 
let us show the solution h N,s of the discretized equation (I3.13P to the solution h 5 for 
the regularized equation (13.111) . when the initial data is smooth enough. The goal is the 
following: 

Theorem 3.16. Assume h e C£°(D). Then, the sequence of solutions {h N ' s (t)} for 
the discretized equation (13.131) with initial datum converges as N — >■ 00 to the unique 
solution h s (t) of (13.111) with initial data ho in the following sense: 

lim \\h N - s {t)-h s {t)\\ m{D) . = Q 

holds for every t > 0. 

Proof. To simplify notations, we omit the parameter S when no confusion arises. 

We shall at first show that we can take a subsequence {N'} such that h N converges to 
the solution of (13.111) . We arbitrarily choose / G C£°(D) such that 



h (6)d6= / f(9)d0 

D JD 

and define / by 

f N (6) = N d [ f(0) d9, 9 e B(x/N, 1/N), x e Z° 



JB(x/N,l/N) 

We introduce the polilinear interpolation used in [2], that is, h is defined by follows: 



h N (t,6)= 

ae{0,l} d 



J(a l {iV^} + (l-a l )(l-{iV^}) 



i=i 



rUm+z), (3 ,8) 
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where [•] and {•} denote the integral and the fractional parts, respectively. We also define 
fc by the similar mannar. We then have 



and 



SUp SUp \\h (t) ||,f/i(_D) < oo 
N 0<t<T 



SUp SUp ||fc (t)||_ffi(D) < oo 
N 0<t<T 



by Proposition ^ . 71 and Corollary 13. 131 Using Proposition [3TZI Corollary 13. 9\ Corollary 13. 131 
and the bounds stated above, we can choose a subsequence {N'} such that 

h N ' ~f N ^9 strongly in C([0, T\,H), 

h N ^h weakly in L 2 ( [0,T\,H% (£>)), 

k N -»■ fc weakly in L 2 ([0, T]xD), 

k N ^ k weakly in L 2 ([0, T], H X (D)) 

as A/ 7 — 7- oo for some g, /i, fc, fc. Letting h = g + f, we can easily see that h = h and 
fc = fc. Furthermore, in this setting, U/i^H jv converges to as A" — > oo. Applying 

the argument in Step 3 of Proposition 1.2 in [5], we obtain that the limit h is the solution 
of (13. lip with initial data Tiq. Furthermore, the uniqueness for (13. lip implies that the 
sequence {h N ;N > 1} itself converges to h strongly in C([0, T], which shows 

the conclusion. □ 



4. Identification of equilibrium states 

In this section, let us study the structure of the equilibrium states for the dynamics 
on (Z d )* corresponding to (12. ip . We will focus our attention to the relationship between 
stationarity and Gibbs property, since we have already known that the family of extremal 
canonical Gibbs measures coinsides with the family of extremal grandcanonical Gibbs 
measures introduced by [5], see [6] for details. 

4.1. Notations. In order to characterize the equlibirum states, we shall prepare several 
notations precisely. Note that we will follow the same manner as in [5] and [TJ. 

Let (7j d )* be the set of all directed bonds b = (x,y), x,y 6 Z d , \x — y\ — 1 in Z d . We 
write Xb = x and = y for b = (x, y). We denote the bond (e^, 0) by e, again if it doesn't 
cause any confusion. For every subset A of Z d , we denote the set of all directed bonds 
included A and touching A by A* and A*, respectively. That is, 

A* := {b E (Z d )*; Xb E A and y b E A}, 

A* := {6 G (Z d )*; x b G A or y b G A}. 

For <p = {4>(x); x G Z d } G M. zd , the gradient V is defined by 

V0(6) :=<t>{x)-<j>{y), b=(x,y)e(Z d )*. 

Now, let X be the family of all gradient fields r] G which satisfy the plaquette 

condition (2.1) in [5], i.e., X = {rj = V</>; G IF*}. Let L 2 . be the set of all rj G R {zd) * 
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such that 

\ v \ 2 r := \v(b)\ 2 e~ 2r ^ < oo. 

We denote X r = X flL 2 . equipped with the norm | • | r . We introduce the dynamics r\ t G X 
governed by the SDEs 

dr) t (b) = -VU.(r) t )(b) dt + V2dVw? (b), b G (Z d )*, (4.1) 

where {wf d (x); x G Z d } is the family of Gaussian processes with mean zero and covariance 
structure 

E[wf{x)wf{y)} = —A Z d(x, y)t A s, x,y G Z d ,t,s> 0. 

Since the coefficients are Lipschitz continuous in X r , this equation has the unique strong 
solution in X r for every r > 0. Note that rjt '■= V</>t defined from the solution <p t of the 
SDE (12. ip on Djy satisfies (14.11) for b G D* N and boundary conditions t]t(b) = Vijj N {b) for 
b G (Z d )* \ D* N , when replacing wf d by w t . 

Let V(X) be the set of all probability measures on X and let V2{X) be those /i G V{X) 
satisfying i? M [|?7(6)| 2 ] < oo for each b G (Z d )*. The measure /i G V^{X) is sometimes 
called tempered. Let Q be the family of translation invariant, tempered Gibbs measures 
/i G V2(X) introduced by [5], and Q cxt be the family of /i G Q with ergodicity under spatial 
shifts. In the case with strict convexity of V, properties of Gibbs measures are studied 
quite well, see [S] and [2]. 

We define the differential operator d x for x G Z d acting on Ci ocb (M. z ) by 

* = -^-. 

<90(x) 

We then have that the generator of (12. ip can be expressed by 

^ = -4^ 9 x (A w a)(x) + 2 ^ (A^t/(V0))(x)9 x . 

We define the differential operator <9 X for x £ Z d acting on Cf oc (X) by 

^ dr,{b) 



We define the differential operator J£ N acting on Cy oc (X) by 

We then have that JzfJv is the generator of r\ t = V(f) t , where 4>t is the solution of (12. ip . 
We also define the differential operator Jzf Z d acting on Cf oc (X) by 

2^ = -4 d x (Ad.)(x) + 2j2(^U(rj))(x)d x . 

x&1 d x&Z d 

To make notations keep simple, we simply denote Jf Z d by Jz? if it does not cause any 
confusion. We also define 

-2a = -4^^(A A a)(x) +2^(A A f/(7 7 ))(x),9 a , 

xeA x£A 
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4.2. Relationship between stationary measures and Gibbs measures. 

Theorem 4.1. Let the probability measure n on X be translation invariant and tempered, 
that is, 

sup E^[r]{b) 2 } < oo 



If fi satisfies 

[ 3?F(r])dn = 
Jx 

for every F G C 2 (X) with compact support, then fi is a canonical Gibbs measure introduced 
in [6]. 

The proof of Theorem 14.11 is similar to [3] , which is based on jl] . 
We at first introduce $a : K — > R by 

$ A ( M )^(l + (Anf)- m , 

where 

a = [ {l + u 2 )- m du. 
it 

For A n := [-n, n] d n 1 d we define $ A : X A * -> R by 

*iw = n $ a(0"'°(x)), 

where r?,a is the height variable satisfying V<f) v,a = r\ and ^'"(O) = a. Note that (j) v,a is 
uniquely determined by r\ and a. We also define Pn(v) by 

pj(»?) = / ^i(i7-0/*w- 

Let = $ A (£ - rj). Since ^(-^) G C 2 C (*), we have 



Multiplying F(£) G Cf oc (X) whose support is in A*, and integrating in £ by the uniform 
measure on Xa* , we obtain 

| F(0^(-,0(»7)A*WdeA. = 0. (4.2) 

Applying Lemma 2.2 in [3] and the relationship 

d^ x (V- Vib) <9$ A <9\l/ A (V</> V0 
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for x G 7L d by the symmetricity of $ A , the left hand side is calculated as follows: 

+ f/fN) £ ^ (-A) (s, y)C/ y fa) ^ J° i/a. W 
=: ii + I 2 , 

where ^A n ,p is the measure on M A ™ defined by 

^a„ iP (#)=p(v(o)) n 

xgA„ 

with a probability density p on ffi. Performing the integration-by-parts for I x , we have 
We shall first calculate I\. Performing integration-by-parts in ip, we have 

- // E(- A )( g ^) F ( v ^ ^ o r ) y ^ (o)) n #o%(^)- (4.3) 

Noting that integrands of l\ and the first term in the right hand side of ( 14. 3 p are function 
of V^, each integral does not depend on the choice of p and therefore the second term 
does not also. Taking a sequence p n such that p' n — > as n — > oo, we conclude that the 
second term must be zero. 
Let us choose F as 

<4 - 4) 

with some bounded smooth function / : K — > K and 

g n fa) = exp (-H A (rj)) , 77 G A^*, 
where Z n is the normalizing constant. Noting 

^n(V-) f dVfaV-) 

d^{x) ~J cty(x) KV} 

and putting 

r x (ri) = 

we have 



xga„ ysA 

+ EE (- A )(^^) / /' KW)) ^ A (V^(WW„, P (#). (4.5) 



zeA„ j/eA 
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Next, we shall compute l^- Performing the integration-by-parts in ip(x) again, we have 



E E(- A )(^) //?!^(^)-^^ 

c&A n y&A n JJ ^ > 

E(- A )^^) JJ^^{u^ 



x&A n y£A 



J2(-^y) II d -^^K^)plw>^ P {d^) 



xeA„ yeA n " " ^(x) 

=: R x (n, f) + R%(n, f) + R x (n, /) + R x (n, f). (4.6) 



Summarizing (14. 3p and f)4.6p . we obtain 



= i? A (n, /) + f) + i? 3 A (n, /) (4.7) 



if we take F as in (14.41) . Note that F x (n, f) is finite when f(x) = logx, by using \A(x, y)\ < 
2d and Lemma 2.3 of [3J. We denote F x (n, f) with f(x) = logx simply by F x (n). 

From now on, we shall show that terms Rf(n, f) in the right hand side can be controled 
by F x (n). 

Lemma 4.2. Assume that the function f satisfies < uf'(u) < 1 for every u > 0. We 
then have bounds for R x (n, f), R x (n, f) and R^(n, f) in (14. 6 j) as follows: 

\R x (nJ)\<KC(n) 1/2 F x (n) 1/2 (4.8) 
|i#(ra,/)| < lfA- 1 F A (fi) 1 / 2 (4.9) 
|i? A (n,/)| < i^(n) 1/2 F A (n) 1/2 (4.10) 

wrf/i a constant K > independent in n and \, where C x (n) is defined by 
C ( n ) = E E / 

xeA n+ i 6g(Z d )*\A*:x,,=a; 

c%(£,n,/*) = / V'(»7(6))A*(^I^X.)(0. 
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Proof. We at first obtain 



/ (£) EE (- A > <■*• v) fi 9r L £L = -f a (") 



dtb(x) dib(y) 



by assumption on /. Note that (—A) is nonnegative definite, indeed, we have 

J](-A)(x,y)0(x)0(y)>O 



for every m > 1 and G M zd such that = on A^. Furthermore, we have the 

Schwarz inequality of the following form: 



1/2 / \ 1/2 



< 



£ £(-A)(x,y)0(x)0(y) £ ^(-A)(x,i/)^(x)^(i/) ) (4.11) 



for every (p,ip E M zd with = ^(x) = on for some m > 1. Using the above, we 
obtain 

for some constant C > 0, which shows ( 14. 8p . We can also obtain (14. 9p and ( l4.1Up by the 
similar argument to the above. Note that (2.12) and (2.13) in [3] is required for the proof 
of (dU). □ 

Let us continue the calculation for F A (n, /). Summarizing (I4.7p - (l4.10p . we get 

F\n, f) < ^F\n) + K'{C{n) + A" V) 

with a constant K' > 0. Using Fatou's lemma, we conclude 

^F\n) < K'(C(n) + A~V). (4.12) 

We shall next give a lower bound for the left hand side of ( I4.12p . For £ £ N, let us take 
A C A„ by 

A= |J At(x), 

where At(x) = Ai + x and 

K,e = {xe ((2£ + 3)Z) d ; A t (x) C A n _!>. 
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Because boxes Ag(x) appearing above are disjoint, we get 

_ I f v V ( - 2 



>-\ j E E E (£g - £g) 

where I A is the entropy production rate defined by 

J A (/i) := sup < / dp,] u G C\ J^r-measurable, u > 1 > 

for a finite A C Z d and /t on or <^vr with n large enough. Applying (I4.12p and taking 
the limit A — > oo, we have 

i £ J A ^(/i)<4K'C(n) 

by the lower semicontinuity of the entropy production rate. Since I Al( ^ x \jjL n ) does not 
depend on x by the translation invariance of \i and C(n) = 0(n d ~ l ), we get 

I Ki (v) < CFrC 1 (4.13) 

with a constant C > indpenent of n. Taking the limit n —> oo, we finally conclude that 

I At (fi) = (4.14) 

for every i G N. Repeating the same argument as in the proof of Lemma 5.2 in [6], we 
obtain that \i is a canonical Gibbs measure introduced in [6]. 

5. Proof of Theorem 12.11 



In this section, we shall give the proof of our main result, Theorem 12.11 Before that, 
we shall prepare several bounds, which play key role in the proof. 

5.1. A priori bounds for stochastic processes. In this section, the goal is to prove 
following proposition: 

Proposition 5.1. There exist constants C\^C>2 independent of N such that 

E\\h N (t)\\\ N + N- d E f J2 (V0f(*)) 2 cfe< ^E]]^ (0)\\\n + C 2 (l+t) 

holds for every t > 0. 
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Proof. Let us use the function / and ip N defined in Section 2.1. We define ifj N by 
where (4>q) is defined by 

r > = j\r*- 1 J2^(y). 

y&D N 



We denote the macroscopic height variable associated with the microscopic height variable 
ip N by f N , that is, 

f N (9) = N-^WlBwwwie), 9 E R d . 

Calculating \\h N (t) — f N \\_iN by Ito's formula, we obtain 

d\\h N (t) - f N \\\ N = —2N~ d (V0f (b) - V^(6)) V{Vtf{b)) dt 

+ 2N~ d i-^y 1 (0f - ^) (*) dw t (y) 

x£D N 

+ 2N- d \D N \dt. 

Therefore, integrating in t and taking the expectation, we get 
E\\h N (T) - f N \\\ N = E\\h N (0) - f N \\\ N 

-2E f N- d V<f?{b)V{V<tf{b))dt 

+ 2E [ N~ d V^ N (b)V'(V<f)?(b))dt 

+ 2N~ d \D N \T. (5.1) 
Applying (13.151) to the third term in the right hand side, we get 



2N~ d ^ N (b)V'(V(f>^(b)) 



1/2 



< 



2c/|(^)l \N- d Yl (nV0f(6))) 2 



beDt 



1/2 



<2c f \(^)\c + lN- d ^(V0f(6)) 2 
<2 lN ' d J2 (V^(6)) 2 + 8 7 - 1 c 2 / |(^)|^ 
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for every 7 > 0. Plugging the above with 7 = c_/2 into (15. ip . we finally obtain 
E\\h N (T)-f N \\\ N <E\\h N (0)-f N \\\ N 

-c.E f N~ d V (V^(b)) 2 dt 

+ 2N- d \D N \T + 16cZ 1 c 2 + c 2 f \(<f ) $)\ 2 T. 
Noting \ (<$) | 2 < \\h N (0) W 2 ^^, we obtain the conclusion. □ 

5.2. Coupled local equilibria. In this subsection, we shall introduce the coupled mea- 
sure and identify its limit point, as in [5] and [7]. Let us denote the (discrete) gradient of 
h N > 5 introduced by (IBTTSj) by uf> a , that is, 

uf' s (x) = V N h N ' 5 (s, x/N), x G Z d , 

and the law of V0f on X by /if . Using these notations, we introduce the coupled measure 
p N (dr]du) on X x R d by 

p N ' 5 (d V du) = t" 1 !^!" 1 I it*s r-^d^s^idu) ds, (5.2) 

-=: — 

where r x is the spatial shift by x G 7L d . We note that the sequence {p N } is tight as the 
probability measures on X x M. d since we have Proposition 15.11 and Proposition 13.71 

Proposition 5.2. For every limit point p of {p N ' s }, there exists a probability measure 
\{dudv) on M. d x M. d such that the relationship 

p(dr]du) = / fj, v (dr))X(dudv) 

holds with ergodic Gibbs measures {// u ; u G M d } introduced by [5]. 

Proof. To keep notation simple, let us omit the parameter 5 when no confusion arises. 
For Lp G Cl(M. d ), we define the signed measure p N (dr], <p) on X by 

p N (dr],ip) = / p N (drj du)<f(u). 

It is sufficient to show that every limit point p(dr), ip) of {p N (dr], <p)} is translation invariant 
and satisfies 

/ &F(r])p(dr],ip) = (5.3) 
for every F G C 2 (X) with a compact support, see Theorem 4.1 of [5]. 
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We shall at first show the limit point p(drj, p) is translation invariant. For F G Cjj(X) 
and e6Z d such that lei = 1, we have 



F(v)p N (drj, tp) — I F(r))p»(;tp)oT e (dr)) 



x 
< 



N i 



.X 



r^jvrl^Hooll^lloo / J2 \u N (s,x/N)-u N (s,x/N -e/N)\ ds 

Jo = 



oo || P || oo 

( \D N n (D N + e) C + \(D N + e) n D N 
Since the first term vanishes as N — > oo by Proposition I3.15[ we obtain 



lim 

N^oo 



F( V )p N (d V ,p)- I F(r l )p N (;p)o Te (dr l ) 



x 



X 



0. 



which shows that the limit p(dr], tp) is translation invariant. 

Let us show (I5.3|) . Fix F G C^{X) with compact support, and take L G N such that 
supp(F) C A* L , where Al = [— L, L] d fl Z d . We then obtain for N large enough 

/ J?F{ V )p N (d V ,p) = r 1 \D N \- 1 [ V p(u^{x))E»° OT *[tfF{rj)]ds 
Jx Jo -=— 

= t~ 1 \D N \- 1 f <P(u?(x))E^[J? DN (FoT^)( V )}ds 

+ t~ 1 \D N \- 1 f ^(x))E^lJ? DN (FoT^)( V )}ds 
Jo 

+ t~ 1 \D N \- 1 f vK(x))E^^[^ DN F(r))]ds 



where the set sd^ l is defined by 



?n,l = D N n I |J (A 2L + z) J 

\zddD N J 



Since supp(F o r x x ) C D N if x G £^n,l, we obtain 

J2 ^(«f(x))^[(For- 1 )(V0f)] 



= t- 1 | J D JV |- 1 AT- 4 



t- 1 


D N \ 


-l N -4 f 






Jo 


r 1 


D N 


-l N -4 f 






Jo 



t 



J2 ^(«f(x))^[(For- 1 )(V0 o Ar )] 
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by Ito's formula. Since F e C%(X), we obtain 



(Ffl < diV" 4 + C 2 N~ 4 f 

Jo 



ds 



with some constants C%,C2 > independent of N. Combining the above with Proposi- 
tion we have 

lim \F?\ = 0. 

For F 2 , we have 



<r l y\UD N \^ f E E^W^For-^Wds 

Jo ^- , .,,« , „ 



<C,\D N \- 1 \D N \ 



+ C 4 \D N \- 1 ! £ £ i^[|V(Vtf(&)|]<fo 

leB^xi/ff.i be(Z d )*:z,,£A 2Z ,+x 

— • r 2,l ^ r 2,2 

with some constants C3, C4 > independent of N. We can easily see that 

lim F£ = 0. 

For i 7 ^, applying Schwarz's inequality, we obtain 

F£ 2 < CsN- x l % \D N \- x E f V |V0f(6)| 2 ds + C % N x ' 2 \D N \- x \dD N \ 

Jo , -p— 



with some constants C 5 ,C e > independent of iV, which indicates 

lim F 2 N 2 = 

N— 5>oo ' 

from Proposition 15.11 We therefore conclude 



lim Frf = 0. 

iV->oc 



Since we obtain 



F 3 N < C 7 N^ 2 \D N \~ l E f Y \V<Ps(b)\ 2 ds + CsN^lDifl-^dDNl 

with some constants C7, Cs > independent of iV by a similar argument to that for F 2 N , 
we obtain 



lim F? = 0. 



Summarizing the above, the identity (15. 3p is concluded. 



□ 
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5.3. Derivation of the macroscopic equation under (13.211) . Let us first prove The- 
orem [2TT] under the following assumption: h N (0) is given by h N (0) = , where Jiq is 
defined by (13TT1) with h G C^(D). Since we have 

\\h N (t) - h(t)\\m { D)* < \\h N (t) - h N > 5 (t)\\ m(Dr 

+ fh N ' 5 (t)-h s (t)\\ m{Dr 

+ \\h s (t)-h(t)\\ HHD) *, 

it is sufficient for our goal to show 

\imsup\imsupE\\h N (t) -h N > 5 (t)\\ 2 = 0, (5.4) 

by using 03. ip . Using Ito's formula, we obtain 
E\\h N (t)-h N >\t)\\\ 



-tN~ 



N 

d „ 
i=l JXxM. d 

+ u i V i a s (u) - l}p N ' s (drjdu) +tN~ d \dD N \, 

where p N ' s is the coupled measure introduced by (15. 2p . Applyin the same argument as in 
the Section 6 of [5] with Propositions 13. 141 and IBTT] we conclude conclude (15. 4p . 

5.4. Derivation of the macroscopic equation in general cases. Let us remove the 
assumption imposed at Section 15.31 and complete the proof of Theorem 12.11 For this aim, 
we prepare the following lemma: 

Lemma 5.3. Let <f) t and <f) t be the solution of (12.11) with common Gaussian processes 
{wt{x); x G Djy} and let h N and h be the macroscopic height variables corresponding to 
<f) t and <pt, respectively. Then, for every t > and N > 1 

E\\h N (t) - h N (t)\\\ N < E\\h N (0) - h N (0)\\\ N 

holds. 

Noting that <p t and <pt satisfy the same boundary condition, the quite same argument 
as in the proof of Proposition 15.11 can be applicable. We therefore omit the proof. 
We shall approximate h G L 2 (D) by h e G Cq(-D) in the sense of 

lim \\h e - h \\ H i {D y = 0. 

Let <f>\ be the solution of (12. ip with common Gaussian processes {w t {x)\ x G D N } as <p t 
and with initial data 0q which are defined by 



0o 0) 



o, xez d \ D N . 
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Letting h N,e be the macroscopic height variable corresponding to <p e t and h E be the solution 
of fll.3p with the initial data h e , we then have 

E\\h N (t) - h{t)f Hl{Dr < AE\\h N {t) - h^{t)f m{Dy + AE\\h^(t) - h%t)f m{D y 

+ s\mt) - h(t)\\ 2 HHDr 

< 4E\\h N (t) - h N m\-i,N + m\h N '%t) - h%t)f m{Dy 

+ 8\mt)-h(t)\\ 2 m{D) , 

by (13.11) . Here, applying the result proved in the previous subsection, Lemma 15.31 and 
Proposition 13.51 we conclude Theorem 12.11 in general settings. 
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